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ABSTRACT 

We present high spectral resolution (R w 108, 000) Stokes V polarimetry of the Classical T Tauri 
stars (CTTSs) GQ Lup and TW Hya obtained with the polarimetric upgrade to the HARPS spec- 
trometer on the ESO 3.6 m telescope. We present data on both photospheric lines and emission lines, 
j—j . concentrating our discussion on the polarization properties of the He I emission lines at 5876 A and 

6678 A. The He I lines in these CTTSs contain both narrow emission cores, believed to come from near 
the accretion shock region on these stars, and broad emission components which may come from either 
a wind or the large scale magnetospheric accretion flow. We detect strong polarization in the narrow 
^vq ' component of the two He I emission lines in both stars. We observe a maximum implied field strength 

J> , of 6.05 ± 0.24 kG in the 5876 A line of GQ Lup, making it the star with the highest field strength 

£Nj ■ measured in this line for a CTTS. We find field strengths in the two He I lines that are consistent 

with each other, in contrast to what has been reported in the literature on at least one star. We do 
not detect any polarization in the broad component of the He I lines on these stars, strengthening the 
conclusion that they form over a substantially different volume relative the formation region of the 
narrow component of the He I lines. 
Subject headings: accretion, accretion disks — line: profiles — stars: atmospheres — stars: formation 
(T) ■ — stars: magnetic fields — stars: pre-main-sequence — 

1. INTRODUCTION Basri 2003, Jayawardhana et al. 2003). Independent of 

T Tauri stars (TTSs) are young ( < 10 Myr), low mass the constraint imposed for defining a CTTS, this 

« , i < n r tvt j. ,i , i i ,i jr moniker has become synonymous with a low mass pre- 

( £ 2.5 Mm stars that have only recently emerged from . , , , , . ,. , J ,. , . , c 

rrt , \_ • iii i ij j. i- i • ii . mam star that is actively accreting material from a cir- 

. v*_ tneir natal molecular cloud cores to become optically vis- , nn , T , , . r , n . n 

r^, , n i j- cumstellar disk. Indeed, the vast majority of stars which 

ible These young, low mass stare are generally subdi- fit ^ criteria for CTTSg snQw gome w Qf additional 

vided into categories such as classical and weak libs. • , r • D n ■ ,• cn , ,. 

™ , . .. j- i -i mmo /nmmm • • n evidence e.g. inverse P-Cygm line profile shapes, opti- 

The designation of a classical TTS (CTTS) was originally , .,. L , b , , x . f . i- j- ( ,. , 

, j , , ,. 1 j- ,1 • cal veiling (sec below), infrared excess indicative of disk 

based on a purely observational distinction: the equiva- accretion 

lent width of the Ha emission line. Classical TTSs are T , . L , . . c ■ 

n^mct i-ii tt -i •1,1 rjr / T t \ 10 It 18 now generally accepted that accretion of circum- 

1 1 Ss which have an Ha equivalent width W eo [Ha) > 10 iiij-ij.-i4.iv, c c ^^^0 • 

■ H 9 v ' stellar disk material onto the surface of a C 1 1 S is con- 

A as distinguished from the weak line TTSs (WTTSs) de- trolled by a strong stellar magnetic field (e.g. see review 

fined by Herbig and Bell (1988); however, Bertout (1989) by Bouvier et al. 2007). These magnetospheric accretion 

suggests that a break point value of 5 A is more appropri- models assert that strong stellar magnetic fields truncate 

ate. More recently, investigators have tied the definition the inner disk, typically near the corotation radius, and 

to the shape (width) of the Ha line profile (e.g. White & channel the accreting disk material onto the stellar sur- 

1 — ttut— : t 7 r 7T -. , face, most often at high stellar latitude (Camenzind 1990; 

cm.i@nce.edu, wc2@nce.edu T ... '. . nnl ri n- 9 u 11 mno 01, I 

vaientMstsci.edu — Komgl 1991; Collier Cameron & Campbell 1993; Shu et 

ieffers@astro.physik.uni-goettingen.de al. 1994; Paatz & Camenzind 1996; Long, Romanova, 

piskunovSStysast.uu.se, Ulcg.Kochukhov@fysast.uu.se, vitaly.makagani^®groa.)k.coffl2(qrj5.^tcr^ Snetohvdrodvnamic 
snik@strw.leidcnuniv.nl, kcllcr@strw.lcidcnuniv.nl, rodenhuis@strw.lcidcnuniv.nl 
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simulations find that outflows launched from near the 
region in which the stellar field interacts with the sur- 
rounding accretion disk can also spin the star down to 
observed rotation rates (e.g. Ferreira 2008; Romanova et 
al. 2009), though some recent work challenges the notion 
that these outflows can actually balance the spin-up ac- 
cretion torques in CTTS systems (e.g. Zanni & Ferreira 
2009). 

Despite the successes of the magnetospheric accretion 
model, open issues remain. Most current theoretical 
models assume the stellar field is a magnetic dipole with 
the magnetic axis aligned with the rotation axis. How- 
ever, recent spectropolarimetric measurements show that 
the fields on TTSs are probably not dipolar (Johns-Krull 
et al. 1999a; Valenti & Johns-Krull 2004; Daou, Johns- 
Krull, & Valenti 2006; Yang, Johns-Krull, & Valenti 
2007; Donati et al. 2007, 2008, 2010a, Hussain et al. 
2009). Few studies of accretion onto CTTSs have taken 
into account non-dipole field geometries. The earliest 
of these by Johns-Krull and Gafford (2002) found that 
abandoning the dipole assumption reconciled observed 
trends in the data with model predictions; however, this 
study did not consider the torque balance on the star and 
whether an equilibrium rotation rate could actually be 
achieved. Johns-Krull and Gafford (2002) argued that 
while the field on the stellar surface may be quite com- 
plex, the dipole component of the field should dominate 
at distance from the star where the interaction with the 
disk is taking place. This assumption appears to gener- 
ally hold true in several recent studies (e.g. Johns-Krull 
& Gafford 2002; Mohanty k Shu 2008; Gregory et al. 
2008; Long et al. 2008; Romanova et al. 2011; Cauley et 
al. 2012). However, the complex nature of the field near 
the surface has significant implications for the size of ac- 
cretion hot spots, making them smaller than would be 
predicted by pure dipole models (Mohanty & Shu 2008; 
Gregory et al. 2008; Long et al. 2008); and also has im- 
portant consequences for disk truncation radii and the 
computation of the torque balance on the star by the 
disk (Gregory et al. 2008; Long et al. 2008; Romanova 
et al. 2011).' 

Two approaches are generally used to measure mag- 
netic fields on low mass stars, both utilizing the Zeeman 
effect. Magnetic fields can be measured from the broad- 
ening of magnetically sensitive lines observed in inten- 
sity spectra (e.g. Johns-Krull 2007; Yang et al. 2008). 
This technique is primarily sensitive to the magnetic field 
modulus, the unsigned value of the field weighted by the 
intensity distribution of the light emitted over the vis- 
ible surface of the star. While this method does not 
suffer from flux cancellation due to regions of opposite 
polarity appearing on the star, it does require that all 
non-magnetic broadening mechanisms be accurately ac- 
counted for in the observed spectra. As a result, this 
technique is primarily sensitive to relatively strong fields. 
Observations of circular polarization in Stokes V spectra 
can be much more sensitive to weak fields on the sur- 
face of stars; however, the Stokes V signature is sensitive 
only to the line of sight component of the magnetic field 
and the signal can be reduced significantly due to flux 
cancellation when opposite field polarities are observed 
simultaneously on the stellar surface. Dopplcr shifts due 
to stellar rotation can reduce the degree of flux cancel- 
lation that results, permitting Stokes V signatures to be 



present even when the net flux weighted line of sight 
field integrated over the stellar surface (the net longitu- 
dinal magnetic field, B z ) is zero. Observations of time 
series of Stokes V spectra can be used to track changes 
in the amount of net field visible on the star as it ro- 
tates, ultimately allowing the large scale field of the star 
to be mapped using various tomographic imaging tech- 
niques (e.g. Donati et al. 2007 and references therein; 
Kochukhov et al. 2004 and references therein). 

In addition to potentially mapping the surface field 
on accreting young stars, information can be obtained 
on the large scale field controlling the interaction of the 
star with its disk and the accretion flow by measuring 
time series of Stokes V profiles in emission lines formed 
in the accretion flow and shock. The first accretion line 
for which circular polarization was detected is the He I 
line at 5876 A (Johns-Krull et al. 1999a), and time se- 
ries of the polarization variations in this line have been 
used to estimate the latitude of accretion spots on sev- 
eral CTTSs (e.g. Valenti & Johns-Krull 2004; Yang et 
al. 2007; Donati et al. 2008, 2010b, 2011a, 2011b). This 
line is observed in most CTTSs and is often found to be 
composed of two components: a narrow core component 
and a broad component extending out to several hun- 
dred km s _1 (e.g. Edwards et al. 1994; Batalha et al. 
1996; Alencar & Basri 2000). Based on the similarity in 
shape between the observed line profiles of some CTTSs 
and model profiles calculated in the context of magne- 
tospheric accretion, Hartmann et al. (1994) suggested 
that the He I 5876 A line (broad and narrow compo- 
nents) might form throughout the accretion flow, with 
the narrow component primarily coming from the lower 
velocity regions near the disk truncation point. Bcris- 
tain et al. (2001) instead argue that the narrow core 
of the He I line arises in decelerating post-shock gas on 
the stellar surface at the base of the accretion footpoints. 
Beristain et al. (2001) argue that the broad component 
observed in many CTTSs has a dual origin in the mag- 
netospheric flow and in a high velocity wind in the most 
strongly accreting stars. 

The strong, ordered fields observed in the narrow com- 
ponent of this line component (~ 2.5 kG, e.g. Johns- 
Krull et al. 1999a) argue for a formation region close to 
the stellar surface instead of several stellar radii above 
the star where the field interacts with the disk. The He I 
5876 A arises from a triplet state and is composed of 
several closely spaced lines. The He I 6678 A line arises 
from the analogous singlet state, and is observed in many 
CTTSs as well where it displays both broad and narrow 
components (see Beristain et al. 2001). Based on the 
strong similarity in their kinematic properties and the 
measured triplet-singlet flux ratio, Beristain et al. (2001) 
conclude that the narrow component of both He I lines 
forms in the post-shock gas. On the other hand, this pic- 
ture is complicated by the observation of Donati et al. 
(2008) that the 6678 A line consistently shows a longitu- 
dinal magnetic field strength approximately twice that of 
the 5876 A line in the CTTS BP Tau whose He I lines are 
dominated by a narrow component (Edwards et al. 1994, 
Batalha et al. 1996, Beristain et al. 2001). This is a sur- 
prising observation since models of accretion shocks on 
CTTSs find that the thickness of the post-shock region is 
typically 10 5 -10 6 cm (Calvet & Gullbring 1998; Lamzin 
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1998) which is a small fraction ( 10~ 5 ) of a stellar ra- 
dius. It would be surprising if the stellar magnetic field 
strength varied so strongly with depth, suggesting then 
that perhaps the two He I lines do not trace the same 
regions on the stellar surface. 

To better clarify the magnetic held properties of accre- 
tion related lines, more spectropolarimetric observations 
of CTTSs, including those with substantial broad compo- 
nents to their He I lines, are needed. Here, we report new 
observations of two CTTSs (GQ Lup and TW Hya) us- 
ing the newly commissioned polarimeter operating with 
the HARPS spectrograph on the ESO 3.6 m telescope 
at La Silla. TW Hya is a K7 CTTS and a member of 
the loose TW Hydrae association (Kastner et al. 1997). 
The Hipparcos parallax for TW Hya implies a distance of 
56 ± 7 pc (Wichmann et al. 1998), making it the closest 
CTTS to the Earth. Based on its placement in the HR 
diagram, the age of TW Hya is estimated to be 10 Myr 
(Webb et al. 1999, Donati et al. 2011b). Setiawan et al. 
(2008) claimed the detection of a ~ 10 Mj up planet in 
a very close orbit around this CTTS, making TW Hya 
an important benchmark constraining the timescale of 
planet formation. Huelamo et al. (2008) instead suggest 
that the observed radial velocity variations which signal 
the presence of the planet are in fact caused by large 
starspots on the surface of TW Hya. As a result, there 
is great interest in knowing as much about this star as 
possible. In addition, TW Hya is still accreting material 
from its circumstellar disk and is observed at a low in- 
clination (i ~ 18°, Alencar & Batalha 2002), making it 
an excellent object for studying magnetically controlled 
accretion onto young stars. The magnetic properties of 
TW Hya have been investigated a number of times pre- 
viously (Yang et al. 2005, 2007; Donati et al. 2011b). 
GQ Lup is also a K7 CTTS, and has also recently come 
under a great deal of scrutiny as the result of a claimed 
planetary mass companion. Neuhauser et al. (2005) dis- 
covered an infrared companion at a separation of <~ 0."7 
(corresponding to ~ 100 AU at a distance of 150 pc). 
Based on their infrared photometry and K band spec- 
tra, Neuhauser et al. (2005) constrained the mass of GQ 
Lup B to be between 1 — 42 Mj up , placing it possibly 
in the planet regime. More recent spectroscopic studies 
have favored the upper end of this range, suggesting the 
companion is more likely a brown dwarf (Mugrauer & 
Neuhauser 2005; Guenther et al. 2005; McElwain et al. 
2007; Seifahrt et al. 2007; Marois et al. 2007; Neuhauser 
et al. 2008). The formation of such an object presents 
challenges to theories of companion formation in a disk, 
and has sparked continued study of this system to bet- 
ter pin down the properties of both of its members. GQ 
Lup is known to show clear signs of variable accretion 
(Batalha et al. 2001), making it a good target to study 
the role of magnetic fields in the accretion process. To 
our knowledge, no studies of the magnetic properties of 
GQ Lup exist to date. In §2 we describe our observa- 
tions and data reduction. The magnetic field analysis 
and results are described in §3, and in §4 we discuss the 
implications of our findings. 

2. OBSERVATIONS AND DATA REDUCTION 

All spectra reported here were obtained at the ESO 3.6 
m telescope on La Silla using the newly comissioned po- 
larimeter, HARPSpol (Snik et al. 2008, 2011; Piskunov 



et al. 2011), mounted in front of the fibers feeding 
the HARPS spectrometer (Mayor et al. 2003). While 
HARPSpol can also record Stokes Q and U spectra, for 
the observations reported here, only Stokes V spectra 
were obtained. As mentioned above, linear polarization 
in both the lines and the continuum can result from scat- 
tering off a circumstellar disk (e.g. Vink et al. 2005); 
however, the action of a disk does not typically produce 
circular polarization in either the lines or the continuum. 
Here, we will focus only on Stokes V in the lines mea- 
sured relative the continuum which is assumed to not be 
circularly polarized. With this instrumental setup, each 
exposure simultaneously records the right and left cir- 
cularly polarized components of the R = 108,000 spec- 
trum. These two components of the echelle spectrum are 
interleaved, such that two copies of each echelle order are 
present on the two 2148 x 4096 CCD arrays (one for the 
blue portion of the spectrum and one for the red). The 
two polarized components of each order are separated by 
^16 pixels in the cross dispersion direction on the array, 
while each spectral trace is ~ 3.5 pixels wide (FWHM) 
in the cross dispersion direction. Each observation of a 
star reported here actually consists of 4 separate observa- 
tions of the star, with the angle of the quarter waveplate 
in the polarimeter advanced by 90° between the expo- 
sures. The result of this is to interchange the sense of 
circular polarization in the two beams. This gives sub- 
stantial redundancy in the analysis which allows us to 
remove most potential sources of spurious polarization 
due to uncalibrated transmission and gain differences in 
the two beams. As described below, we use the "ratio" 
method to combine the spectra from these interchanged 
beams in order to form Stokes / and V spectra that are 
largely free of these potential spurious signals (e.g. Do- 
nati et al. 1997; Bagnulo et al. 2009). All spectra were 
obtained on the nights 29 April 2010 through 2 May 2010, 
with one night (1 May) lost due to weather. TableQ]gives 
a complete table of the stellar observations reported here. 
Included in the Table are continuum signal-to-noise es- 
timates near the two He I emission lines studied here as 
well as the emission equivalent widths of these two lines. 
Also reported is the veiling found near the He I 6678 A 
line as discussed below. Along with spectra of GQ Lup 
and TW Hya, a spectrum of the weakly accreting TTS 
V2129 Oph was also obtained and is used in the analysis 
of the He I lines on the other stars. In addition to stellar 
spectra, standard calibration observations were obtained 
including bias frames, spectra of a Thorium Argon lamp 
for wavelength calibration, and spectra of an incandes- 
cent lamp for the purpose of flat fielding. The calibration 
spectra were obtained with the polarimeter in front of the 
fibers. 

All spectra were reduced with the REDUCE package 
of IDL echelle reduction routines (Piskunov & Valenti 
2002) which builds on the data reduction procedures de- 
scribed by Valenti (1994) and Hinkle et al. (2000). The 
reduction procedure is quite standard and includes bias 
subtraction, flat fielding by a normalized flat spectrum, 
scattered light subtraction, and optimal extraction of the 
spectrum. The blaze function of the echelle spectrometer 
is removed to first order by dividing the observed stel- 
lar spectra by an extracted spectrum of the flat lamp. 
Final continuum normalization was accomplished by fit- 
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ting a 2nd order polynomial to the blaze corrected spec- 
tra in the regions around the lines of interest for this 
study. Special care was taken to apply a consistent con- 
tinuum normalization procedure to the spectra extracted 
from all four sub-exposures. Occasional small difference 
in normalization of the two orthogonal spectra are com- 
pensated by using the "ratio" method (e.g. Bagnulo et 
al. 2009, and below) to combine the right and left circu- 
larly polarized components. The wavelength solution for 
each polarization component was determined by fitting 
a two-dimensional polynomial to nX as function of pixel 
and order number, n, for approximately 1000 extracted 
thorium lines observed from the internal lamp assembly. 
The resolution as determined by the median FWHM of 
these thorium lines was R = 107, 660. 

As mentioned above, each subexposure obtained of a 
given star contains both the right and left circularly po- 
larized component of the spectrum. In order to get a hnal 
measurement of the mean longitudinal magnetic field, 
B z , these individual measurements of the two circular 
polarization components must be combined in some way. 
We used the "ratio" method (e.g. Bagnulo et al. 2009, 
Donati et al. 1997) to combined the right and left circu- 
larly polarized components of the spectra form the Stokes 
V spectrum as well as a null spectrum, with each being 
renormalized to the continuum intensity. We also added 
all the components together to form the Stokes / spec- 
trum. With Stokes V and / determined, the continuum 
normalized right-hand circularly polarized (RCP) com- 
ponent of the spectrum is then R = I + V and the con- 
tinuum normalized left-hand circularly polarized (LCP) 
component of the spectrum is L = I — V. Computing 
these from / and V in this way ensures both circular po- 
larization states have been normalized to the same con- 
tinuum. 

3. ANALYSIS 

3.1. He I Line Equivalent Widths 

Table[T] gives the equivalent width of the two He I lines 
studied here for all our target stars. As mentioned before, 
previous investigators have noted that these lines often 
appear to have two distinct components (e.g. Batalha 
et al. 1996; Beristain et al. 2001): a narrow component 
(NC) and broad component (BC). It is thought that the 
two components may form in different physical regions 
of the accretion flow onto CTTSs (Beristain et al. 2001) 
and their polarization properties also appear to be dif- 
ferent with the narrow compnent showing significantly 
stronger polarization (Daou et al. 2006; Donati et al. 
2011b). We therefore report the equivalent width of the 
NC and the BC separately for the two He I lines, the 
sum giving the total line equivalent width. Decomposing 
the lines in this manner requires some assumptions to be 
made about how to separate the two components. Since 
the NC often appears asymmetric (e.g. Figure 1) with 
a very steep blue edge and shallower red edge, Gaussian 
fitting to the lines requires particular choices to be made 
on just how to do the analysis. For example, Batalha et 
al. (1996) define (by eye) a region outside the NC and fit 
a single Caussian to the resulting BC and subtract it off 
in order to measure the NC equivalent width. Another 
procedure is to fit the entire line with multiple Gaussians 
and use the resulting fit parameters to estimate compo- 



nent properties (e.g. Alencar & Basri 2000). The result- 
ing equivalent with of the various features then depedns 
at some level on how one chooses to do the analysis. This 
is illustrated in Figure 1. The top panel shows the He I 
5876 A line of TW Hya from the first night. The smooth 
solid curve shows a line profile fit employing 3 Gaussian 
components. The dash-dot line shows a fit using only two 
Gaussian components. There is a clear difference in the 
two fits (the 3 Gaussian fit uses two Gaussians to fit the 
NC which is not really Gaussian as mentioned above). 

The bottom panel of Figure 1 zooms in on the line 
to show the recovered BC profiles. The BC from the 3 
Gaussian fit is shown in the smooth solid line and that 
from the two Gaussian fit is shown in the dash-dot line. 
Also shown is a BC fit (dash-triple dot line) following 
Batalha et al. (1996) where a single Gaussian is used to 
fit the region outside the NC. Finally, the solid straight 
line connecting the two large squares shows a by eye es- 
timate of the point on both the blue and red side of the 
line (as seen in Stokes /) where the NC and BC join with 
a linear interpolation between these points to define the 
separation of the NC and the BC which can be used to 
separately determine their equivalent widths. This then 
gives 4 different ways to estimate the equivalent width 
of the BC (and also the NC). The two extremes for the 
BC equivalent width are the single Gaussian fit (1.193 
A) following Batalha et al. and that (1.089 A) from the 
two Gaussian fit, corresponding to a difference of 9%. 
Clearly, none of the Gaussian fits exactly follow the red 
side of the BC, so it is impossible to predict just what this 
component does under the NC. Given this uncertainty 
and the fact that using the linear interpolation between 
the blue and red sides of the apparent boundary between 
the NC and BC gives equivalent width values between the 
two extremes, we choose to use this method to separate 
both components and measuure the equivalent widths re- 
ported in Table [1] We note that the BC equivalent width 
for the profile shown in Figure 1 computed this way dif- 
fers from that resulting from the 3 Gaussian fit by only 
4.9%. We therefore estimate that the systematic unc- 
etainty resulting from the choice of just how to separate 
the two components likely leads to a 5% uncertainty in 
the reported equivalent widths which is not included in 
the Table. In most cases, the boundary between the BC 
and the NC is clear and repeated measurements with 
slightly different choices yield results with a difference 
less than la for the quoted uncertainties. There are a 
few cases where the boundary between the NC and the 
BC, or the BC and the continuum, are less clear and we 
repeated the measurements with a larger distinction in 
our choices of these points. These are noted Table [Hand 
we use our different measurement trials to estimate the 
the equivalent width uncertainty for these profiles. For 
the other measurements, the uncertainties are computed 
by propagating the uncertainties in the observed spectra. 

3.2. The Photospheric Mean Longitudinal Field 

For each of the T Tauri stars, we measured the photo- 
spheric B z using approximately 40 magnetically sensitive 
absorption lines (Table 2), which form primarily over the 
portions of the stellar surface that are at photospheric 
temperatures. These lines may have relatively little con- 
tribution from the cool spots that are likely present on 
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these stars. Due to the wavelength dependence of the 
Zeeman effect and the fact that the signal-to-noise ra- 
tio achieved in the observations of these late-type stars 
is considerably higher in the red regions of spectrum, 
we focus the analysis here only on the spectra from the 
red CCD of HARPSpol. Lines for the analysis were se- 
lected by visual inspection of all the orders on the red 
CCD recorded with HARPS. Lines were deemed good 
for the analysis if they appeared relatively strong (cen- 
tral depth £ 0.15) in the observed spectrum (though 
most were considerably stronger), appeared free of blend- 
ing by other photospheric lines, and were not contami- 
nated by telluric absorption. Lines passing these criteria 
were then checked in the Vienna Atomic Line Database 
(VALD, Kupka et al. 1999, 2000) and if they are present 
and have a value for the effective Lande g-factor, the 
line was used in the analysis. In a few cases, the VALD 
data indicated that an apparently good line is actually 
a very close blend of two lines. In this case, we used 
the line but estimated a new effective Lande (/-factor by 
calculating the weighted mean of the effective Lande g- 
factors of the lines in the blend. The weights used are 
the central depth of each component line as predicted 
by VALD for the atmospheric parameters typical of K7 
TTS (T eff = 4000 K, log.g = 3.5). The initial line list 
was constructed using a visual examination of the spec- 
trum of GQ Lup obtained on 29 April 2010. For the other 
TTSs some lines were affected by blending with telluric 
absorption or by strong cosmic ray hits (as is also the 
case for later observations of GQ Lup). In these cases, 
the lines were not included in the determination of the 
photospheric B z values. Lines so affected are noted in 
Tabic 2. 

Once the line list was determined, the mean longitu- 
dinal magnetic field, B z , can be estimated by measuring 
the wavelength shift of each line, AA = \r — Al, where 
Xji is the wavelength of the line observed in the RCP 
component of the spectrum and Al is the wavelength 
measured in the LCP component of the spectrum (Bab- 
cock 1962). The shift of the line observed in the two 
polarization states is related to B z by 

AA = 2 A £ 9 X 2 g cB B z = 9.34 x IQ~ 7 \ 2 g cB B z mA 

47T7TI e Cr 

(1) 

where g e s is the effective Lande g-factor of the transi- 
tion, B z is the strength of the mean longitudinal mag- 
netic field in kilogauss, and A is the wavelength of the 
transition in Angstroms (Babcock 1962; see also Mathys 
1989, 1991). In order to measure the wavelength shift, 
AA, we measured the wavelength of each line of interest 
in the two circular polarization components using the so- 
called center of gravity method (e.g. Mathys 1989, 1991; 
Plachinda & Tarasova 1999). This method for estimating 
B z is mathematically equivalent to estimating B z from 
the first order moment of the Stokes V profile, assuming 
the same integration limits are used in the two methods 
(e.g. Borra & Vaughn 1977; Mathys 1989; Landi Deglln- 
nocenti 2004). Determining B z requires knowledge of 
the effective Lande g- factor, g e g, for the transition. The 
weights for individual 7r and a components of a given 
spectral line which go into the definition of g e s assume an 
optically thin line, so equation (1) is only approximately 
true in the case of moderately strong (saturated) photo- 



spheric lines. As a result, using equation (1) is strictly 
valid only for weak lines, but give good results for real 
spectral lines (Mathys 1991). For each measurement of 
the center of gravity of a spectral line, we used the locally 
measured signal-to-noise in the spectrum to estimate the 
uncertainty in the spectrum and then used standard er- 
ror propagation to find the uncertainty in the line shift 
bewteen the two polarization states and the implied B z . 
Our final estimate of B z is a weighted mean of the in- 
dividual line estimates, and these means are reported in 
Table 3. 

We repeated the measurements of the line wavelengths 
in the two polarization components and the resulting fi- 
nal value of B z a number of times, making slightly differ- 
ent choices on integration limits for the center of gravity 
estimate of the wavelength of each line. However, for 
each trial we always used the same integration limits for 
a given line when analyzing the RCP and LCP compo- 
nents of the spectra. In each case, we achieved consistent 
results within our quoted uncertainties. Generally, we di- 
vided our trials in two groups. In the first case, we choose 
integration limits very close to where the lines appear to 
reach the local continuum. This was primarily done as 
an effort to exclude any potential weak line blends that 
might appear as a small distortion in the line wings. In 
the second group, we choose integration limits clearly out 
in the local continuum, but which in some cases likely in- 
cluded some weak line blends. In many cases, choosing 
the wider limits produced higher values of B Zl though 
in some cases the measured field went down. On aver- 
age, the wider bins resulted in fields stronger by ~ 35%. 
Choosing the wider limits does generally result in some- 
what larger uncertainty estimates, typically by a factor 
of 1.8. Again, the two groups of results are consistent 
within these uncertainties (the difference typically being 
~ 1.5c). In Table 3 we quote the values for the wider 
integration limits with their correspondingly greater un- 
certainties. 

Examining the photospheric fields and their uncertain- 
ties as reported in Table 3, apparently significant fields 
are found on both TW Hya and GQ Lup each night they 
were observed. However, the value of B z found on V2129 
Oph is less than la and does not represent a real detec- 
tion. As a test of our measurement techniques and in 
order to gain confidence in our uncertainty estimates, 
two different null tests were performed on the observa- 
tions of each target. Each test should return a measured 
value of B z = 0, thus serving to test how accurately we 
can recover this value and whether the uncertainties are 
realistically estimated. As described above, the spectra 
from the 4 subexposures were also combined in such a 
way as to produce a null Stokes V profile (e.g. Donati 
et al. 1997; Bagnulo et al. 2009), which can be used to 
calculate null RCP and LCP spectra. We analyzed these 
null spectra in the same as described above. This has 
the advantage of using exactly the same lines as used to 
measure B z , with exactly the same wavelength limits for 
computing the shift of each line, but with data combined 
in a different way than in the real measurement. The sec- 
ond null test we employed used 15 strong telluric lines 
from the atmospheric B band of oxygen between 6883 - 
6910 A. Since these lines should show no significant po- 
larization of their own, we can analyze them in exactly 
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the same fashion as we do the stellar lines from which 
we derive B z (that is we take A A = \r — Xl as done for 
the stellar measurements) combining the data from the 
sub-exposures in exactly the same way as done for the 
stellar measurements. In order to translate the measured 
shifts to a value of B z we assign a Lande ^-factor to each 
telluric line equal to the mean value of the photosphcric 
lines used for the given observation. The weights are the 
uncertainty on the value of B z derived from each of the 
photospheric lines. These null test field values are also 
reported in Table 3. 

3.3. B z in the Accretion Shock Emission 
As described above, Johns-Krull et al. (1999a) dis- 
covered that the He I 5876 A emission line can be cir- 
cularly polarized in spectra of CTTS, implying coherent 
magnetic fields at the footpoints of accretion columns. 
They measured B z = 2.46 ± 0.12 kG for BP Tau. Since 
this original discovery, polarization in the He I 5876 A 
emission line has now been reported for a number of 
CTTSs by a number of investigators (Valenti & Johns- 
Krull 2004; Symington et al. 2005; Smirnov et al. 2004; 
Yang et al. 2007; Donati et al. 2007, 2008, 2010b; Chen 
& Johns-Krull 2011). Polarization has since been de- 
tected in other emission lines (notably He I 6678 A and 
the Ca II IRT lines) thought to be associated with accre- 
tion shock emission as well (Yang et al. 2007; Donati et 
al. 2007, 2008, 2010b; Chen & Johns-Krull 2011). The 
spectra obtained here contain both the He I 5876 A and 
6678 A lines, so we analyze them with a focus on trying 
to see how well the fields derived from each line agree 
with one another as this could provide clues to the loca- 
tion and geometry of the accretion shocks on the star. As 
mentioned earlier, the 5876 A line of He I is composed 
of several components (6) which are closely spaced in 
wavelength, and as a result is subject to the Paschen- 
Back effect (e.g. Yang et al. 2007; Asensio Ramos et 
al. 2008). Therefore, the exact splitting pattern of the 
lines can vary considerably, depending on the strength 
of the magnetic field. However, since most of the level 
crossing and merging has occurred by the time the local 
field strength reaches 2 kG (e.g. Asensio Ramos et al. 
2008), the He I 5876 A line should be in or very close to 
the complete Paschen-Back limit given the field strengths 
we recover below. As a result, we set g e ff = 1.0 for this 
line as done in Yang et al. (2007). We can then test 
the validity of this assumption once we have our field 
measurements . 

Figure 2 shows the right- and left-circularly polarized 
components of the He I 5876 A emission lines of GQ Lup 
and TW Hya as observed on 29 April 2010. Also shown 
in the figure are the Stokes V profiles of the lines. Nearby 
photospheric absorption lines are also seen in each star. 
Since the polarization signal in the photosphere is quite 
weak due to the low value of B z present there (Table 
3), these individual photosphcric absorption lines do not 
show obvious polarization. They do serve to show that 
the two polarization components are well aligned in wave- 
length, so that the obvious shift of the He I emission line 
between the RCP and LCP components indicates a very 
strong field in the line formation region. 

In order to measure the value of B z in the line for- 
mation region of the He I line we again use the center 



of gravity technique to measure the wavelength shift of 
the line as observed in the two polarization components 
(RCP and LCP). Our measurements of B z and its uncer- 
tainty for both He I emission lines are reported in Table 
3 for each star on each night. When measuring the field 
in the He I formation region, care must be given when se- 
lecting the wavelength limits for determing the center of 
gravity wavelength of the line and also in deciding how to 
separate the narrow (NC) and broad (BC) components 
of the line (§3.1). Looking at the He I 5876 A line of TW 
Hya in Figures 1 and 2, polarization is clearly seen in the 
NC of the line, but is not obviously apparent in the BC 
extending off to the red side of the line. As a result, we 
focussed in on the NC of the line when measuring the 
value of B z in the 5876 A line of He I. 

In addition to the specific wavelength region chosen, 
care must also be taken when defining the local contin- 
uum to be used when measuring the center of gravity 
wavelength for the line. The reason for this is that the 
center of gravity technique (as well as the first moment 
of Stokes V) is an intensity weighted mean wavelength, 
where the intensity used is that above the continuum in 
the case of an emission line. In the case of the He I 5876 
A emission line, the NC of the line sits on top of a BC in 
many cases as discussed above. In this case, significantly 
different results are obtained if the stellar continuum is 
used compared to what is obtained if a somewhat higher 
continuum defined by the BC is used. We proceed in 
an effort to isolate the emission from the NC and mea- 
sure B z in this component. Interpreting the NC of the 
He I emission as an excess line emitted from a distinct 
region that adds its emission to that from both the stel- 
lar continuum and the BC of the He I emission, these 
additional sources should be subtracted off when mea- 
suring the center of gravity wavelength of the NC of the 
emission line. In §3.1 we described several methods of 
separating the NC and BC when measuring their equiv- 
alent width, showing that each method is subject to cer- 
tain biases, but that the resulting systematic differences 
were small (~ 5% for the method we settle on). We used 
the same methods to remove the BC from the line (each 
BC is shown in the bottom of Figure 1) and measured 
the resulting B z . For the profile shown in Figure 1, us- 
ing the single Gaussian fitted to the region outside the 
NC gave the largest magnitude field (-2.57 ± 0.11 kG), 
while the two Gaussian fit and the linear interpolation 
both gave —2.38 ± 0.10 kG, and the 3 Gaussian fit gave 
—2.40 ± 0.11 kG. All 4 methods give results consistent 
to within 2a and three of the 4 differ by 0.02 kG or less. 
As in the case for the equivalent widths above, we again 
adopt the linear interpolation under the NC as the way of 
removing the BC and use the resulting NC to determine 
the B z values given in Table 3. 

As mentioned earlier, there have been reports in the 
literature that the He I 6678 A line shows substantially 
stronger polarization than the 5876 A indicative of a 
stronger local field in the line formation region of this 
more optically thin line. We looked for evidence of this 
effect in our stars by analyzing the the 6678 A emission 
line in all 3 of the T Tauri stars oberved here in a way 
similar to how the the 5876 A line is analyzed. How- 
ever, it became immediately apparent that additional 
care needs to be taken when analyzing the 6678 A line. 
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Figure 3 illustrates the situation for GQ Lup with the 
spectrum obtained on the first night of our observing 
run. Shown in the figure is the right and left circularly 
polarized components of the spectrum, along with the 
Stokes V profile. Unlike the 5876 A line where the right 
and left circularly polarized components of the emission 
line have essentially the same shape, only shifted, the 
6678 A line of GQ Lup has a different shape in the two 
circular polarization components. On the other hand, 
the Stokes V profile for the two He I lines looks quite 
similar. The reason this is the case is the presence of an 
Fe I photospheric absorption line at the same wavelength 
of the 6678 A line which can be seen in the spectrum of 
V2129 Oph (taken on night 2) which is also shown in 
Figure 4. In V2129 Oph, the He I 6678 A line shows 
no obvious emission, but it turns out there is some weak 
emission in this line in V2129 Oph that partially fills in 
the photospheric absorption line and produces a clear 
Stokes V signature which is illustrated in Figure 4 and 
discussed below. 

As described above, the formation region for the nar- 
row components of these He I lines is thought to be at 
the base of the accretion columns where material from 
the disk is raining down onto the star. In addition to 
producing some emission lines, these accretion footpoints 
are believed to be the source of the optical continuum 
veiling seen in most CTTSs (e.g. Basri & Batalha 1990; 
Hartigan et al. 1991; Valenti et al. 1993; Gullbring et al. 
1998). This extra line and continuum emission region ef- 
fectively blocks some small portion of the stellar surface 
and the light it emits, adding its own emission on top 
of the stellar spectrum coming from the non-accreting 
regions of the star (e.g. Cavet & Gullbring 1998). The 
spectrum of the excess can then be studied by subtract- 
ing an appropritely scaled (veiled) stellar template spec- 
trum from the observed CTTS spectrum (e.g. Hartigan 
et al. 1995; Gullbring et al. 1998; Stempels & Piskunov 
2003). The spectrum of the excess is sometimes stud- 
ied by subtracting off the spectrum of a veiled WTTS 
or main sequence star of the same spectral type, while 
other studies use a synthetic spectrum computed from a 
model stellar atmosphere. Here, we use a model stellar 
atmosphere since we did not observe a suitable WTTS 
having no excess continuum or line emission of its own. 

The top panel of Figure 3 shows a portion of the spec- 
trum of V2129 Oph in the neighborhood of the He I 6678 
A line with our best fit synthetic spectrum including a 
modest amount of veiling (r — 0.17, that is, a continuum 
excess equal to 0.17 of the local stellar continuum) added 
in. This value for the veiling is slightly higher than previ- 
ously reported values for V2129 Oph: r = 0.0±0.1 (Basri 
& Batalha 1990); r = 0.075 ±0.020 (Donati et al. 2011a, 
though this is only a relative veiling value and represents 
a lower limit to the true value). For our purposes, the 
exact value of the veiling is not important. The primary 
goal with such a fit is to predict the strength of the Fe I 
absorption line at 6678 A. To do so, we synthesize the 
spectral region between 6660 - 6680 A as shown in Fig- 
ure 4. The two strongest features in this region are the 
Fe I feature at 6678 A and another Fe I feature at 6663 A. 
Both of these features seen in Figure 4 are actualy blends 
of two Fe I lines. However, since all the components of 
both features are from the same element in the same 



ionization stage, it should be possible to "fit" the 6663 
A feature and then predict the strength of the 6678 A 
line. We use the spectrum synthesis code SYNTHMAG 
(Piskunov 1999). The line data for this spectral region 
is taken from VALD (Krupka et al. 1999, 2000). As is 
often the case, the initial predicted line strengths do not 
match up well with observations of the TTSs or the Sun, 
so we tuned the oscillator strengths (and in a few cases 
the Van der Waals broadening constants) of the strong 
lines until a good match with the solar atlas (Kurucz et 
al. 1984) was obtained. In tuning the line parameters 
to the solar spectrum, we use the model atmosphere and 
associated parameters (vsim, macroturbulence, etc.) for 
the T e f f = 5731 K scaled Kurucz model from Valenti and 
Piskunov (1996). Once the line data is set, we then syn- 
thesize this same spectral region for each of our CTTSs 
using NextGen model atmospheres (Allard & Hauschildt 
1995). To do so, we must select an effective tempera- 
ture, gravity, and metallicity for each star. We take the 
gravity as logg = 3.5 and metallicity as [M/H] = 0.0 for 
all stars; however, we pick effective temperatures from 
the standard NextGen grid that are as close as possible 
to that for each star. Donati et al. (2007) adopt an 
effective temperature of 4500 ± 200 K for V2129 Oph, 
and we adopt the NextGen model with T e ff = 4600 K 
for this star. Yang et al. (2005) find T eff = 4126 ± 24 
K for TW Hya, and we adopt the NextGen model with 
T e ff — 4200 K. GQ Lup has the same spectral type as 
TW Hya (K7), so we use the same NextGen model for 
this star as well. In order to do the final spectrum syn- 
thesis and fitting to the observed profile, we must select 
the micro and macroturbulence as well as the wsini rota- 
tional velocity. Following Johns-Krull et al. (1999b) we 
take 2.0 km s _1 for the macroturbulent velocity and set 
the microturbulence to 1.0 km s _1 . In order to account 
for rotation, we set wsim = 4.0 km s _1 for TW Hya (Do- 
nati et al. 2011b), we set vswi = 14.5 km s" 1 for V2129 
Oph (Donati et al. 2007), and we set vsini = 6.8 km s _1 
for GQ Lup (Guenther et al. 2005). 

Using the spectrum synthesis described above, we 
match the strength of the Fe I feature at 6663 A and 
subtract the model spectrum from the right and left cir- 
cularly polarized components of the observed spectrum, 
add a pseudo-continuum of 1.0 back in, and then follow 
the procedure described above for the He I 5876 A line in 
order to measure the field in the formation region of the 
6678 A He I line. The only free parameter used to fit the 
6663 A Fe I is the value of the continuum veiling. These 
veiling values and the B z values for the 6678 A emission 
line are reported in Table 1 . In general our veiling values 
agree well with previously published determinations con- 
sidering this quantity is often quite variable in CTTSs. 
In addition to previous veiling measurements for V2129 
Oph discussed above, veiling on TW Hya in similar spec- 
tral regions has been shown to vary, reaching as high as 
0.80 in the study of Alencar and Batalha (2002) and as 
high as <~ 0.92 in the work of Donati et al. (2011b). In 
the case of GQ Lup, Weise et al. (2010) found a veiling 
of 0.5 near the Li I line at 6707 A, while here we find 
values ranging from 0.3 to 0.6. Again though, the exact 
value of the veiling is of secondary importance: account- 
ing for the photospheric absorption which is coincident 
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with the He I 6678 A emission makes a substantial dif- 
ference in the recovered field strengths for this line. To 
gauge this effect, we repeated the analysis of the 6678 A 
line without making any correction for the photospheric 
absorption present. Generally, the fields we measure in 
this case are a factor of two larger than those reported 
in Table 3. 

4. DISCUSSION 

Examining both the null Stokes V field values and the 
field values obtained from the analysis of the telluric 
lines shows that each of these is quite low and generally 
equal to zero (as they should be) within the measured 
uncertainties. In the case of the null spectrum, the value 
reaches a significance of 3.3(7 in one case, but most values 
are between 1 — 2o\ For the telluric lines, all of the mea- 
sured B z values are equal to zero (as expected) to within 
la or less. The uncertainties from the telluric tests are 
also generally lower than the for the photospheric lines 
(both actual measurement and null test). This is primar- 
ily due to the telluric lines being very sharp and strong, 
allowing for very accurate measurements of any potential 
wavelength shift. We conclude from these tests that our 
uncertainties in B z are generally well characterized and 
that our detections of polarization in the photospheric 
lines and resulting measurements of B z on GQ Lup and 
TW Hya are real. 

TW Hya has been studied with spectropolarimetry by 
Yang ct al. (2007) and Donati et al. (2011b). On one 
of six nights, Yang et al. (2007) measured polarization 
in a dozen photospheric lines of TW Hya correspond- 
ing to a longitudinal magnetic field of B z = 149 ± 33 
G, while finding no significant polarization (though with 
larger uncertainties) on the other five nights. Donati et 
al. (2011b) observed TW Hya for a total of 20 nights 
and used their least squares deconvolution (LSD) anal- 
ysis (Donati et al. 1997) to measure photospheric fields 
ranging from B z = 380 — 700 G with uncertainties typi- 
cally of 15 G. Donati et al. suggest that long term tem- 
poral variability may be responsible for the variations in 
B z between the measurements of Yang et al. (taken in 
April 1999) and those of Donati et al. (taken in March 
2008 and March 2010). The results presented here, from 
April and May 2010, agree well with those of Yang et al. 
and less so with those of Donati et al. (2011b). 

Yang ct al. (2007) measured polarization in the He I 
5876 A and 6678 A emission lines of TW Hya, as well as 
in the Ca II 8498 A emission line. Yang et al. report pos- 
sible rotational modulation of the polarization in the He I 
5876 A line, with implied fields varying from ~ —1450 to 
~ —1800 G. Donati et al. (2011b) also measure strong 
polarization in the He I 5876 A line with variations again 
suggestive of rotational modulation. The field strength 
implied in this case ranges from ~ —2000 to <~ —3500 
G. Donati et al. attribute the significant difference be- 
tween their He I measurements and those of Yang et al. 
to perceived differences in the way the fields were mea- 
sured in the two studies. Donati et al. measured the 
field only in the narrow component of the line for simi- 
lar reasons to those cited above in §3.3, and they state 
that Yang et al. (2007) used the entire He I line in their 
field determination; however, this is not correct. Yang et 
al. (2007) also focussed only on the narrow component 



of the He I emission line in their field determinations 
(Yang, private communication) . The field measurements 
reported here in Table 3 for B z in the He I 5876 A emis- 
sion line generally agree well with those of Yang et al. 
and are significantly less than most of the fields reported 
by Donati et al. However, Yang et al. used a cross cor- 
relation technique to measure line shifts and resulting 
field strengths, while here we use the center of gravity 
technique to measure line shifts. 

In order to verify that the difference in measurement 
technique does not introduce a spurious difference in the 
recovered field, we reanalyzed the Yang et al. (2007) 
spectrum of TW Hya obtained on 21 April 1999 (the 
strongest He I field they found) using the same center of 
gravity technique employed above. Briefly, this spectrum 
was obtained with the 2.7 m Harlan J. Smith telescope 
of McDonald Observatory used to feed the Robert G. 
Tull coude echelle spectrometer (Tull et al. 1995). The 
starlight was split into its circularly polarized compo- 
nents using a Zeeman Analyzer (ZA, Vogt et al. 1980) 
placed in front of the spectrometer slit. The ZA con- 
tains a Babinet-Soleil phase compensator used to correct 
for a potential reduction in sensitivity to circular polar- 
ization which can be introduced by the non-azimuthally 
symmetric reflections in the coude mirror train. Such 
bounces can introduce some linear polarization into an 
originally circularly polarized beam and the phase com- 
pensator is used to convert this linear polarization back 
into the original circular polarized signal (see Vogt et 
al. 1980 and Vogt 1978). The exposure of TW Hya for 
that night totalled 4300 s. Yang et al. (2007) measured 
a field of —1806 ± 114 G for the narrow component of 
the He I 5876 A line using a cross correlation analysis, 
while we find a field of —2326 ± 118 G using the center of 
gravity technique, fully consistent with our current mea- 
surements. This measurement of the field is stronger by 
520 ± 164 G than that of Yang et al. (2007), represent- 
ing a 3.2cr difference. The difference is driven both by 
the method of estimating the line shift (center of grav- 
ity versus cross correlation) as well as in the treatment 
of removing the broad component. As a result, it does 
not appear that differences in the measurement technique 
can account for the variations in the field strength in the 
He I emission line formation region found between Do- 
nati et al. (2011b) and this study plus that of Yang et 
al. (2007). It is more likely that intrinsic variability is 
at work in this accretion diagnostic, and indeed, Donati 
et al. (2011b) find significant differences in the He I field 
strength from one rotation phase to the next in TW Hya, 
while the bulk of their field measurements are larger in 
magnitude than either Yang et al. (2007) or those here. 

One of the motivations for this study was to verify 
and expand on the intriguing result that the field mea- 
sured from the He I 6678 A line is significantly different 
(stronger) than that measured in the He I 5876 A line 
(e.g. Donati et al. 2008). In principle, measuring the 
field in different emission line diagnostics could offer a 
means for probing the magnetic field structure through 
the accretion shock on CTTSs. As described in §3.3, 
we discovered that there is a strong photospheric Fe I 
absorption line coincident in wavelength with the He I 
6678 A line that can severely affect the field strength 
measured in this emission line if the photospheric line 
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is not properly accounted for. As described above, we 
attempted to correct for the Fe I photospheric line by 
computing veiled model spectra for each of our observa- 
tions and subtracted the resulting model from the obser- 
vations before measuring the held. Doing so produced 
field measurements in the He I 6678 A line that are very 
consistent with those measured in the 5876 A line of the 
same clement. The field strengths recovered from the 
two lines are the same to within 3<r for all observations. 
If we do not account the Fe I photospheric absorption 
line, we generally recover a field in the 6678 A line that 
is twice as strong as reported in Table 3. In a separate 
test, we corrected for the photospheric absorption line in 
the observations of TW Hya and GQ Lup using a veiled 
version of the observed V2129 Oph spectrum. While this 
is not ideal due to the weak He I 6678 A emission from 
V2129 Oph, the results were the same as those in Ta- 
ble 3 within the uncertainties. As a result, it appears 
that the field in the two He I lines is essentially identi- 
cal, at least in TW Hya and GQ Lup, but that special 
care must be taken when analyzing the 6678 A line. We 
also note from earlier that our analysis of the 5876 A 
line assumed it formed in the complete Paschen-Back 
regime (i.e. g e ff = 1.0). Since shock models suggest 
the field strength in these two lines should essentially be 
the same, and that is indeed what we find, these results 
also suggest that our treatment of the 5876 A line in the 
complete Paschen-Back effect is appropriate, at least for 
the strong magnetic fields recovered here. 

The magnetic field on GQ Lup appears to be quite 
remarkable. The field detected in the photospheric ab- 
sorption lines is fairly typical in magnitude (\B Z \ ~ 200 
G) relative to the strength detected on many CTTSs. 
However, the field (~ 6 kG) detected in the He I 
line formation region, the accretion shock, is substan- 
tially stronger in magnitude than that observed in most 
CTTSs. Symington et al. (2005) reported some quite 
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strong fields in this line from a few CTTSs; however, also 
with substantial error bars. Higher precision measure- 
ments (with uncertainties from ~ 0.1 — 0.4 kG) have peak 
fields measured in the He I 5876 A line of \B Z \ < 3.4 kG 
(Donati et al. 2011) with peak values typically \B Z \ ~ 2 
kG depending on the specific star in question (Johns- 
Krull et al. 1999a; Valenti & Johns-Krull 2004; Yang 
et al. 2005; Donati et al. 2007, 2008, 2010b). Very 
recently, Donati et al. (2012) present polarimetric mea- 
surements of GQ Lup at two different epochs (July 2009, 
June 2011). While the photospheric fields they detect 
do show some variation from one epoch to the other, our 
measurements over about half the rotation period found 
by Donati et al. are consistent with the fields found at 
both epochs. On the thcr hand, Donati et al. find a 
substantial decrease in the large scale field controlling 
the accretion over the two years between their epochs. 
Our data (April 2010) is fully consistent with their ear- 
lier epoch, with our He I B z determinations matching 
the strongest values they observe in July 2009, providing 
additional constraints on the timescale involved in the 
magnetic field change observed by Donati et al. 
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Table 1 

Observing Log 





UT 




S/N h 


S/N b 


He I 5876A 


He I 5876A 


He I 6678A 


He I 6678A 




UT Date 


Time a 


Star 


5876A 


6678A 


NC W eq (A) 


BC W eq (A) 


NC W eq (A) 


BC W eq (A) 




29 April 2010 


1:00 


TW Hya 


59 


56 


1.110 ±0.004 


1.179 ± 0.007 


0.387 ±0.003 


0.269 ± 0.005 


1.00 




4:16 


GQ Lup 


42 


42 


0.604 ± 0.005 


0.137 ± 0.008 


0.129 ±0.003 


0.030 ± 0.005 


0.40 


30 April 2010 


0:56 


TW Hya 


66 


78 


1.282 ± 0.004 


2.037 ± 0.008 


0.363 ± 0.002 


0.479 ± 0.004 


1.40 


3:02 


GQ Lup 


31 


29 


0.530 ± 0.005 


0.235 ± 0.009 


0.175 ± 0.005 




0.65 




7:28 


V2129 Oph d 


65 


68 


0.134 ±0.002 








0.175 


02 May 2010 


3:09 


TW Hya 


90 


82 


1.312 ± 0.003 


1.307 ±0.005 


0.393 ±0.011 c 


0.281 ±0.011 


1.00 




6:35 


GQ Lup 


83 


72 


0.365 ± 0.002 


0.620 ± 0.020 f 


0.037 ± 0.006 


0.110 ±0.007 


0.30 



a This is the midpoint of the 4 X 1800 s exposures that make up each total observation. 

b This is the S/N in the continuum near the respective He I emission lines, calculated from the 

final Stokes / spectrum. 

c This is the veiling in the vicinity of the 6678 A He I line. 

d The entries for V2129 Oph which show no data are due to either there being no clear broad 
component in the case of the 5876 A line, or there being no emission above the continuum in 
the case of the 6678 A line. As discussed later in the text, there is some filling in of a nearby 
photosphcric absorption line by He I emission at 6678 A; however, the overall line remains below 
the continuum level and so we do not record and emission equivalent width here. 
e Due to apparent photosphcric absorption on the blue side of the line, there is some ambiguity 
in how to separate the NC and the BC on this side of the line profile. The reported value and 
larger uncertainty here takes into account repeated measurements (averaging to get the value) 
where more or less of the emission is attributed to the BC or the NC. 

f The line on this night appeared to have quite extended BC wings, making it difficult to establish 
exactly where the line rejoined the continuum. The measurement and uncertainty are formed 
by averaging conservative and more broadly inclusive measurements of the 5876 A line for this 
night. 
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Table 2 

Lines Used for Photosphcric Field Analysis 



Element Wavelength (A) Lande-g e ^ 



V I 


6058.139 


2.14 


Ti I 


6064.626 


1.99 


Fe I 


6173.334 


2.50 


Blend 


6216.355 


1.59 


Fe I 


6219.278 


1.66 


Fe I d 


6232.640 


1.99 


Fe I 


6246.316 


1.58 


V I 


6251.827 


1.57 


Fe I 


n o cr o cr cr A 


0.95 


V I 


6274.648 


1.53 


V I 


6285.149 


1.58 


Or I 


6330.091 


1.83 


Fe I 


6336.823 


2.00 


Fe I 


6393.600 


0.91 


Fe I 


6408.018 


1.01 


Fe I 


6411.647 


1.18 


Fe I 


6421.349 


1.50 


Ca I a 


6439.075 


1.12 


Blend 


6462.629 


0.98 


Ca I 


6471.662 


1.20 


Fe I b 


6475.624 


1.90 


Fe I 


6481.869 


1.50 


Ca I 


6493.780 


0.88 


Fe I 


6498.937 


1.38 


Ca I c 


6499.649 


0.96 


Fe I 


6518.365 


1.15 


V I 


6531.415 


1.57 


Cr I a 


6537.921 


1.71 


Fe I 


6574.227 


1.25 


Ni I 


6586.308 


1.02 


Fe I 


6593.870 


1.13 


Ti I 


6599.105 


0.98 


VI 


6624.838 


1.43 


Ni I 


6643.628 


1.31 


Fe I 


6663.334 


1.53 


Li I 


6707.799 


1.25 


Fe I 


6710.316 


1.69 


Cal 


6717.681 


1.01 


Ti I 


6743.122 


1.01 


Fe I 


6750.149 


1.50 



a This line excluded from all TW Hya analysis due to 
apparent additional blending. 

b This line excluded from analysis of GQ Lup on 30 April 
2010 and from analysis of V2129 Oph due to significant 
cosmic ray hit. 

c This line excluded from analysis of GQ Lup on 02 May 
2010 due to significant cosmic ray hit. 
d This line excluded from analysis of V2129 Oph due to 
significant cosmic ray hit. 
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Table 3 

Magnetic Field Measurements 



UT Date 


Star 


B z (Phot) 
(G) 


B z (Null) 
(G) 


B z (Tel.) 

(G) 


B z (He I 5876A) 
(kG) 


B z (He I 6678A) 
(kG) 


29 April 2010 


TW Hya 


181 ±27 


8 ±27 


-3 ±7 


-2.38 ±0.10 


-2.46 ± 0.13 




GQ Lup 


-195 ±35 


0±35 


-1 ± 11 


6.05 ± 0.24 


6.15 ±0.39 


30 April 2010 


TW Hya 


176 ± 20 


-15 ±20 


2 ± 5 


-1.87 ±0.09 


-2.14 ±0.11 




GQ Lup 


-236 ± 49 


-52 ± 49 


-4 ±14 


5.63 ±0.39 


6.39 ±0.50 




V2129 Oph 


-19 ±32 


-54 ±31 


-8 ±8 


2.87 ±0.39 


3.10 ±0.56 


02 May 2010 


TW Hya 


180 ± 16 


-13 ± 16 


0±6 


-2.25 ±0.09 


-2.53 ±0.08 


GQ Lup 


-81 ± 18 


6 ± 18 


-4 ±6 


5.23 ±0.18 


4.87 ±0.31 
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Figure 1. In each panel, the continuum normalized He I 5876 A line profile of TW Hya from 29 April 2010 is shown in the solid histogram. 
The top panel shows two multi Gaussian fits to the profile, with a two Gaussian fit shown in the dash-dot line and a 3 Gaussian fit shown 
in the smooth solid line. The bottom panel zooms in on the line to show the recovered BC profiles. The BC from the 3 Gaussian fit is 
shown in the smooth solid line and that from the two Gaussian fit is shown in the dash-dot line. The dash-triple dot line shows a single 
Gaussian fit following Batalha ct al. (1996). The solid straight line connecting the two large squares shows a by eye estimate of the point 
on both the blue and red side of the line (as seen in Stokes /) where the NC and BC join with a linear interpolation between these points 
to define the separation of the NC and the BC. 
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Figure 2. In each panel, the continuum normalized right circularly polarized component of the He I 5876 A line profile is shown in the 
upper black curve and the left circularly polarized component is shown in the gray curve. The difference (right - left) of the two, offset by 
0.5, is shown in the lower black curve. This difference is 2x Stokes V with the profiles normalized in this way. 
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Figure 3. The continuum normalized right circularly polarized component of the He I 6678 A line profile of GQ Lup is shown in the 
upper black curve and the left circularly polarized component is shown in the gray curve. The difference (right - left) of the two, offset by 
0.2, is shown in the bottom black curve. This difference is 2x Stokes V with the profiles normalized in this way. The middle black curve 
shows the same region in the Stokes / spectrum of V2129 Oph, revealing the presence of a strong Fe I photosphcric absorption line which 
can alter the inferred field in the He formation region if the Fe line is not properly accounted for. 
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Figure 4. The upper panel shows the observed Stokes I spectrum of V2129 Oph in an expanded region near the He I 6678 A line. The two 
strongest lines are Fe I lines at 6678 A and 6663 A. The smooth gray curve in the upper panel shows a synthetic fit to the spectrum, showing 
that there is excess emission the spectrum of V2129 Oph at 6678 A due to a weak He I emission line superimposed on the photosphcric 
absorption spectrum. In the bottom panel, the upper black curve shows the continuum normalized right circularly polarized component of 
the He I 6678 A line profile of V2129 Oph after subtracting off the synthetic fit and adding f .0 back to the result. The gray curve shows 
the same for the left circularly polarized component of the spectrum. The difference (right - left) of the two, offset by 0.8, is shown in the 
lower black curve. This difference is 2x Stokes V with the profiles normalized in this way. 



